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i — Abstract A simulation of electric current pulses formed by a packet of gamma-quanta moving 
'^hrough an absorptive medium is presented. The electromagnetic fields of the current pulse 
i moving along the straight line with super light velocity are obtained. 
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Introduction 



>in recent years the electromagnetic fields of sources moving with velocities equal or higher than 
*^elocity of light c in vacuum have been discussed in connection with an interest to the problem of 
i— idealized electromagnetic waves [1, 2]. Such a source can be obtained by interaction of a directed 

_beam of a hard electromagnetic radiation (X-rays) with an absorptive rod. It is important to have 

>the length of this rod larger than its diameter. If initial quanta are moving along the axis of this 
^•od, the electron's current pulse whose velocity is equal to c is generated. If the beam is directed at 
^^n angle to the rod's axis, the current pulse has a super light velocity. The velocities of electrons, of 
^ourse, are smaller than c, but the motion of area, where the current's density is non zero, can be 
ly-presented as the motion of some effective charge which is faster than the velocity of light in vacuum. 
3T'he investigation of the spatial structure and of the time dependence of secondary electromagnetic 

fields, generated by such a pulse, requires solution of two problems. Firstly, we have to find the shape 

• the current pulse generated by the beam of primary hard electromagnetic radiation. Secondly, we 
>4^ave to solve the Maxwell's equation for this current's pulse with zero initial conditions. The current 

'■^jyulse in question was earlier studied mainly on a phenomenological basis. The current pulse moving 
^with velocity of light was considered in connection with nuclear explosions [3]. 

• ^4n this paper we calculate the "realistic" shape of current pulse moving with super-light velocity, 
^4king into account the microscopic movements of the electrons. We obtained electromagnetic fields 

cdtvith this current pulse by using the approach presented in [4], based on the method of an incomplete 
separation of variables by V.I. Smirnov [5] and Riemann formula. 



2 The shape of the super-luminal current 

pulse generated by a shot-lived beam of a hard radiation 

The solution of the first problem of finding the shape of current pulse requires to take into account 
the photon absorption's processes, Compton scattering, and secondary coupling effects of "delta 
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electrons" with matter which result in the formation of the ions and the electrons. We found the 
shape of current pulse via the computational experiment using the method of numerical modelling 
described in [6] . This method enables us to study the time dependence of the electrons' distribution 
in the phase space. The code GEANT [7] was used. The generation of the electrons with energy 
lower than 10 keV was not taken into account. 

The scheme and the initial conditions of the computational experiment are as follows. The 
absorbing area is placed in vacuum and is bounded by the cylindrical surface of 0.3 m. in diameter 
and of 50 m. in length L and by the two planes, which are orthogonal to the cylinder axis. Its walls 
are made of hydrocarbon-based material. Their thickness is 0.001 m. The absorptive medium is air 
at the pressure of 1 atm. The symmetry axis of the absorbing area coincides with the axis z of a 
cylindrical coordinate system . 

The origin of the coordinate system is chosen to be localized at the point O, which lies in one of 
the two above mentioned planes. 

The packet of the primary photons (a gamma-ray pulse) starts from a surface formed by rotation 
of a line segment, which is at an angle n/2- 9 to the Z axis . The photons are simultaneously emitted 
at an angle 9 to the direction of the z axis. 

The computational experiment is beginning with the starting time of the gamma-ray pulse. The 
results of computer experiments are as follows: 

1. We obtained the distribution N(p) of the secondary photoelectrons which have passed through 
the plane fixed at z=z . This distribution describes the electron localization in the radial direction 
P- 

2. The width (AR) of the distribution N(p) calculated for the above initial conditions is less 
than 30 centimeters (AR < 30cm). The obtained value meets the necessary requirement of the 
applicability of the model of a linear current (L >> AR), used in the simplified electrodynamic 
calculations. 

The time dependence of the current's pulse passing through the absorptive area at z Q = 30 meters 
is shown in Fig. 1. 

j 




Figure 1: The normalized current pulse in the cross section of the absorbing region at z = 30 m. 
In this plot the value of function in maximum is presented as 1.0 on the axis of ordinates. The 
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pulse length (AT) taken on the ordinates level of 0.1 does not exceed 0.2 ns (or 6 cm in ct - units). 
The obtained value of AT estimates the localization of the current along z - axis. 
The shape of the running pulse can be described by the following function: 

J(z, r) = A • (1 — exp(— ot\((3r — z))) • exp(— a 2 (/3T — z)) (1) 

ai > 0, i=l,2, 

where r = ct is the time variable, (5 — v/c , A is constant ; coefficients ai are obtained by the least 
square method 

Thus, as a result of the above computational experiment the space-time description of the hard 
radiation inducing super light running pulse was obtained. 



3 Calculation of the magnetic component of the electro- 
magnetic wave generated by the current pulse moving 
with super light velocity 

The basic solution of the Maxwell equations that describes waves generated by a current pulse 
propagating along the section of a strait line was given in [2]. In that work the components of both 
the electric and magnetic fields are expressed via a single scalar function v(p,z,r). In the present 
work we represent only magnetic component B. That was made to simplify the subsequent physical 
interpretation of the obtained results. 

In our case the parameters of the hard radiation and absorptive medium, geometry of the ab- 
sorptive area were chosen in accordance with the condition L » AR. The boundary effects in 
the description of a source of electrodynamics problem are not taken into account (the necessary 
requirement L » AT is met). This means that for the calculation of the required electromagnetic 
field (its magnetic component) we may use the simplified model of a linear current. 

Thus the source in the electrodynamics problem (non-zero component of the current's density 
vector) can be taken in the form: 

j z = 5 4^h(PT-z)h(z)h(L-z)J(z,T) r>0 , (2) 
J Z = r<0 

Here S(p) is the Dirac distribution, h(z) is the Heaviside step function. In the case of the axial - 
symmetric source of an electrodynamics problem, the components of the intensity E of electric field 
and the induction B of magnetic field may be expressed via one scalar function v according to the 
following expressions: 

__(Pu_ d 2 u t d 2 u _c?u_ 

p - dpdz z - dri + dz 2 9 ~ drdp ' [ ' 

Now the electrodynamics problem can be expressed in the form: 

( d 2 1 d ( d \ d 2 \ 4vr 

- — « = = <°, (4) 



\dr 2 pdp\ dp J dz 2 J c 

were v = du/dr The non-zero component of the magnetic induction .B^ vector is equal to—dv/dp. 

The solution of the problem (4) in space - time representation is defined by the following expres- 
sion: 
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t'+z+t 



v(p,z,r) = - J dr j dz'h(z') 



t'+z-t 



5(t' - T + Jp 2 + (z- z') 2 ) 
h((3r' - z ')h(L-z r )J(z',r')— y ' 



V 2 + (z- z') 2 

After interchanged the order of an integration, we receive: 

1 f 1 

v(p,z,r) = - / dz'h{z')h{L-z') h 

C z-r \Jp 2 +(z-z') 2 
1 Z+T 1 

+- / dz'h(z')h(L- z') . 1 2 

c J v/p 2 + (z- z') 2 



where 



-z'+T — Z 



h= J dr'J{z', r')h(Pr' - z')5(r' - r + yj ' p 2 + (z - z') 2 ) 



-z'+T + Z 



I 2 = j dr'J{z', r')h{f3r' - z')8(r' - r + yj ' p 2 + (z - z') 2 ) 
o 

For a super light delta - pulse 

Jz = 2^ 6{PT " Z)h{z)k{l ~ z)f{z ' T) 

(where / is differentiable function) expression (6) can be simplified 

v(p,z,r) = -(^1 + ^2) 
c 



Here 



^i(p, z,t) = 

( ( Pz-t + VD \ Pz-t + VD \ 1 - h(y/D + (1 - f3 2 )r - /3(z - L3r)) 

J \ \ P 2 -i )' P 2 -i ) Vd 

x{[h(L)(h(z - t) - h(z + r))] 

x [h(P{VD -zP + r)) - h(-zf3 2 + L(-l + (3 2 ) + (3(V15 + r))] 

- z - r)h(z + r) 

x + + (-1 + + P 2 )t) - h(-z(3 2 + L(-l + /3 2 ) + /3(v^ + r))] 

- r)/i(L - z + r) 

x + VDf3 + (1 + P- P 2 )t) + h(-zp 2 + L(-l + /3 2 ) + /3(v^D + r))]} 



f ( a ( (3z-T-\fD \ Pz-t- VD \ 1 - h(y/D + (1 - (3 2 )t + (3(z - f3r)) 

x/ V \ /5 2 -i ) ' /? 2 -i ) 7B 

x{[h(L)(h(z-r)-h(z + r))\ 

x [h(-(3(VD + z(3- r)) - h(-zf3 2 + L(-l + /3 2 ) + (/? - v 7 ^ + r))} (8) 

- ^ - r)/i(z + r) 

x + VD(3 - (-1 + (3 + (3 2 ))t - h(-z(3 2 + L(-l + (3 2 ) + /3(-v^D + r))] 

- r)/i(L - z + r) 

x + v^D/? + (-1-/3 + (3 2 )t) + h(-z(3 2 + L(-l + (3 2 ) + (3(-Vd + r))]}, 

where D = (z — (3t) 2 — p 2 ((3 2 — 1). This expression may be used also for the nondispersive matter 
(j3 is the ratio of velocity of the source motion to the light velocity in given environment in this case). 

We carried out the numerical calculation of super light source using the expression (6). The time 
dependence of v (p, z, r) is presented on fig. 2 for the seven angles of gamma - quanta. 
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Figure 2: Time dependence (in ct units) of v(p, z, r) normalized to the maximal amplitude of different 
observation angles = 0.3°(1), 9 = 3°(2), 9 = 5.5°,9 = 5.44°(1), 9 = arccos(c/w) = arctg(0.1) w 
5.71°(4), 9 = 6°(5),0 = 9°(6), 9 = 10°(7) 



In the proximity of the angle 9 = arccos(c/t') = 5.71° the function v(p, z, r) sharpens its form. 
Note, that the variation of parameters of numerical experiment causes essential changes of v(p, z, r) 
. In particular, the reduction of value of L results in disappearance of a characteristic maximum in 
the distribution v{9) at 9 = arccos(c/t> ). The dependence of B v component on r at various angles of 
observation 9 (L = 50m, Z = 50000m) is presented in Fig. 3. 

As we see, duration of the two-polar pulses of an electromagnetic field depends on angle 9. 
The reduction of duration and the growth of amplitude of the first half-wave in the vicinity of 
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Figure 3: B v components of magnetic induction at various angles of observation 9 = 0.3°(1), 9 = 3°(2) 
, 9 = arccos(c/t;) = arctg(0.1) « 5.71°, 9 = 6° (A), 9 = 9°(5) 9 = 10°(6) and L = 50m, Z = 50000m 

9 = arccos(c/w ) is observed. The representation B^p, z, if, r) by Fourier integral allowed us to obtain 
the frequency distribution of energy of the electromagnetic field for different values of observation 
angles in a far zone. The frequency spectra of electromagnetic pulses for observation angles : 9 = 0.3°, 
9 = 3° , 9 = 5.5° ,6 = 5.44°, 9 = arccos(c/t;) = arctg(0.1) « 5.71°, 9 = 6°,9 = 9° 9 = 10° in the 
optical range are presented in Fig. 4. 
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Figure 4: The energy spectra of electromagnetic pulse observed at various angles 9 = 6.4°(1), 
9 = 0.3°(2), 9 = 3° , 9 = 5.5° (4) , 9 = 5.8°(5), 9 = arccos(c/w) = arctg(0.1) w 5.71°(6) in the far 
zone ( Z = 50000m) in optical frequency range (3 * 10 14 — 8 * 10 14 ) Hz for v/c= 1.005. 



As we see, radiation at the angle 9 = arccos(c/t>) to the direction of current pulse motion in the 
visible range of a spectrum has maximal value. 



4 Summary 

We have calculated the shape of a the current pulse moving with velocity greater then the velocity 
of light in vacuum for the certain geometry of the hard radiation at fixed values of parameters of 
gamma - quanta pulse moving through a layer of air. 

The analytic representation for the potential v (p, z, r) was obtained for the case of the superlu- 
minal delta pulse. 

The spectra of radiation in the optical range for some different observation angles range for v/c 
equal 1.005 were calculated. We found that the radiation has maximum at 9 = arccos(c/t>). 
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